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Abstract 


The Weibull distribution is a very applicable model for the lifetime data. In this paper, 
we have investigated inference on the parameters of Weibull distribution based on record 
values. We first propose a simple and exact test and a confidence interval for the shape 
parameter. Then, in addition to a generalized confidence interval, a generalized test variable 
is derived for the scale parameter when the shape parameter is unknown. The paper presents 
a simple and exact joint confidence region as well. In all cases, simulation studies show that 
the proposed approaches are more satisfactory and reliable than previous methods. All 
proposed approaches are illustrated using a real example. 
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1 Introduction 

The Weibull distribution is a well-known distribution that is widely used for lifetime models. 

It has numerous varieties of shapes and demonstrates considerable flexibility that enables it 
to have increasing and decreasing failure rates. Therefore, it is used for many applications for 
example in hydrology, industrial engineering, weather forecasting and insurance. The Weibull 
distribution with parameters a and /?, denoted by W (a,/3), has a cumulative distribution 
function (cdf) 



X > 0, a > 0, /3 > 0, 
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and probability density function (pdf) 


/ (x) =, X > 0. 
aP 

The Weibull distribution is a generalization of the exponential distribution and Rayleigh 
distribution. Also, Y = log (X) has the Gumbel distribution with parameters b = ^ and 
a = log (a), when X has a Weibull distribution with parameters a and /3. 

Let Xi,X 2 ,... be an infinite sequence of independent identically distributed random vari¬ 
ables from a same population with the cdf Fg, where 0 is a parameter. An observation Xj will 
be called an upper record value (or simply a record) if its value exceeds that of all previous 
observations. Thus, Xj is a record if Xj > Xi for every i < j. An analogous definition deals 
with lower record values. The record value sequence {Rn} is defined by 


Rn = Xt„, n = 0,1,2,.... 


where is called the record time of nth record and is defined as = min{j : Xj > X'r„_i} 
with To = 1. 

Let Ro,..., Rnhe the first n +1 upper record values from the cdf Fg and the pdf fg. Then, 
the joint distribution of the first n + 1 record values is given by 

n—1 


Jr (r) = fe (rn) J 


1=0 


fe jri) 

-Fg{ny 


ro < ri < • • • < r„, 


( 1 . 1 ) 


whe re r = (rn , 


ri,.. ., Vn) and R = {Rq, Ri,..., Rn) (for more details see 


Arnold et al. 


19981 1. 


Chandlen (j 19521 ) launched a statistical study of the record values, record times and inter 


record times. Record values and the associate d statistics a re of intere st and import a nce in the 
areas of meteorology, sports and economics. 


Ahsanullah 

(1995 

) and 

Arnold et al. 

(19981 


two good references about records and their properties. 

Some papers considered inference on the Weibull distribution based o n record values: 


( 1982 1 discussed some distributional results based on upper record values. 


Dallas 


Balakrishnan and Chan 


(|l994l ) established some simple recurrence relations satisfied by the single and the product mo- 


ments , and derived the BLUE of the scale parameter when the shape parameter is known. 


Chad (j 19981 1 provided a conditional method to derive exact intervals for location and scale 
parameters of location-scale fam ily that can be used to derive exact intervals for the shape 


parameter. 


Wu and Tseng (|2006|) provided some pivotal quantities to test and establish con- 


fidence interva l of th e shape parameter based on the first n -|- 1 observed upper record values. 


Soliman et al. 


(12001)1 ) derived the Bayes estimates based on record values for the parameters 





























with re spect to squared error loss function and LINEX loss func tion. 


Asgharzadeh and Abdi 


compu ted the coefficient of skewness of upper/lower record statistics. 


(j2011bl l proposed joint confidence regions for the parameters. iTeimouri and Guptal (1201^') 


Teimouri and Nadaraiah 


(|2013l i derived exact expressions for constructing bias corrected maximum likelihoo d estimators 


(MLE ’s) of the parameters for the Weibull distribution based on upper records. 


Gouetetal 


(j2012l i obtained the asymptotic properties for the counting process of (5-records among the first 
n observations. 

In this paper, we consider inference about the parameters of Weibull distribution based on 
record values. Eirst, we will propose a simple and exact method for constructing conhdence 
interval and testing the hypotheses about the shape parameter j3. Then using the concepts 
of generalized ]>value and generalized confidenc e interval, a generalized app roach for infer¬ 


ence about the scale parameter a w ill be derivec . 
concept of generalized p-value, and IWeerahandi 


Tsui and Weerahandil (Il989l i introduced the 


(jl99^ i introduced the concept of generalized 


conhdence interval. These 

problems fsee Weerahandi 

approaches have been used successfullv to address several complex 
199511 such as conhdence interval for the common mean of several 

log-normal distributions ( 

Behboodian and 

afari. 

2006 

), conhdence interval for the mean of 

erence about the stress-strength relia- 

Weibull distribution (Krishnamoorthv et ah. 

2001 

), inf 


bility involving two independent Weibull distrib utions 


comparing two dependent generalized variances (j.Tafaril . 


Krishnamoorthv and Linl . 


20ld h and 


2013). 


We also present an exact joint confidence region for the parameters. Our simulation studies 
show that the area of our joint confidence region is smaller than those provided by other existing 
methods. 

The rest of this article is organized as follows: A simple method for inference about shape 
parameter and a generalized approach for inference about the scale parameter are proposed in 
Section [2j Eurthermore, a simulation study is performed and a real example is proposed in this 
Section. We also present a joint conhdence region for the parameters a and /3 in Section [3l 


2 Inference on the parameters 

Suppose Rq, Ri,Rn are the hrst n -|- 1 upper record values from a Weibull distribution with 
parameters a and /3. In this section, we consider inference on the parameters a and (3. Erom 




























































(HI), the joint distribution of these record values can be written as 

< ro < n < • • • < 


i=0 


( 2 . 1 ) 


Therefore, {Rn,Y17=o^^^ sufficient statistic for {a, 13). Moreover, it can be easily shown 

that the MLE’s of the parameters a and f3 are 

n +1 Rr, 




Er=o log V -Ri 


Rn 


a = 


(n + l)/3 


( 2 . 2 ) 


Theorem 2.1. Let Rq, i?i, ..., Rn he the first n + 1 upper record values from a Weibull distri¬ 
bution. Then 


i. 17 = 2/3 X)r=o log ® chi-square distribution with 2n degrees of freedom, 

a. V = 2(-^)^ has a chi-square distribution with 2n + 2 degrees of freedom. 

Hi. U and V are independent. 

Proof, i. Dehne 


Qm — 


Rr, 


From 


and 


Arnold et al. 


Rm—l 


-, m = 1, 2,..., n. 


Therefore, 


(Il998l ) page 20, Qm’s are independent random variables with 

PiQm> q) = q>i, 

2f3mlog{Qm) = 2/3mlog( ) ~ yL- 

n—1 


C/ = 2fi^log{^) = 2f3Y.log( 


i=0 
n—l n 


i=0 




Rti 


Rn—1 Rn—2 

n m—1 


Ri 


Rm—l 


2/3 EE log((5m) = ^2 2;5"^log(Qm), 


(2.3) 


i=0 m=i+l ^ m=l i=0 m=l 

has a chi-square distribution with 2n degrees of freedom, 
ii. Dehne 

X ^ 

^ = (-) , 
a 

where X has a Weibull distribution with parameters a and j3. Then, Y has an exponential dis¬ 
tribution with parameter one. Ther efore, we can conclu de that V has a chi-square distribution 


with 2n -|- 2 degrees of freedom (see 


Arnold et al. 


19981 . page 9). 

















iii. Let (5 be known. Then, it can be concluded from (|2.ip that is a complete sufficient statis¬ 
tic for a. Also, Qm’s in (I2.3P are ancillary statistics. Therefore, Rn and QmS are independent, 
and the proof is completed. □ 

2.1 Inference on the shape parameter 

Here, we consider inference on the shape parameter, /3 from a Weibull distribution based on 
record values, and propose a simple and exact method for constructing a confidence interval 
and testing the one-sided hypotheses 

Ho:/3</3o vs. Hi :/3 > /3o, (2.4) 

and the two-sided hypotheses 

Ho:/3 = /3o vs. ili :^/3o, (2.5) 


where /3o is a specified value. 

Based on Theorem l2.ll U = 2(3 Y17=o (■^) ^ chi-square distribution with 2n degrees 

of freedom. Therefore, a 100 (1 — 7 ) % confidence interval for (5 can be obtained as 

2 2 
^(2n),7/2 A(2n),1-7/2 

2Et„iog(t)’2E".,iog(t) 

where xfk) 7 0'^^ percentile of the chi-square distribution with k degrees of freedom. Also, 

for testing the hypotheses in (|2.4I) and (12.5p . we can define the test statistic 

C/„ = 2/3„glog(^^), 

and the null hypothesis in (|2.4p is rejected at nominal level 7 if 

Uq > X(2n),l-7’ 

and the null hypothesis in (|2.5p is rejected if 

Uq < X(2n),7/2 ^0 > X(2n),l-7/2- 



Wu and Tsend (j2006lf proposed the random variable 

p/3 


Wi/3) = 


(^+i)(nr=o^*)^ 








for inference about the shape parameter, and showed that W{I3) is an increasing function with 
respect to /3. Also, its distribution does not depend on the parameters a and /?. In fact, Vb(/3) 
is distributed as 


= 


E n rj* 

j=0 


(n+i)(nr=o^*)'^ 

where R* is the ith record from the expone ntial distribution with parameter one. However, its 
exact distribution is very complicated, and IWu and Tseng (120061 ) obtained the percentiles of 


W (P) using Monte Carlo simulation. The conhdence limits for /3 are obtained by solving the 
following equations numerically as 


W{(3) = W*_^/„ W{(3) = w;/„ 


where ITJ is the 5th percentile of the distribution of W*. 


(2.7) 


2.2 Inference on the scale parameter 

Here, we consider inference about the scale parameter, a for a Weibull distribution based on 
record values, and propose an approach for constructing a confidence interval and testing the 
one-sided hypotheses 

Hq : a < ao vs. Hi : a > no, (2.8) 

and the two-sided hypotheses 

Hq : a = ao vs. Hi : a ^ ao, (2.9) 


where oq is a specihed value. 

We did not hnd any approach in literature for inference about a based on record values when 


the shape parameter is unknown. Here, we use the concepts of gene ralizec 


alized confidence interval introduced by 


Tsui and Weerahandil (jl989|), and 


jo-value and gener - 


Weerahandi ( 1993 1 


respec tively. In appendix, we briefly review these concepts, and refer readers to 
(I 1995 I I for more details. 

Let 

2C'r _ C'r _ 

r = r„(2)“ 

V itn 


Weerahandi 


( 2 . 10 ) 


where Cr = (^) ’ r*, i = 0,1,... , n is the observed value of Ri, i = 0,1,... , n, 

and U and V are independent random variables that are dehned in Theorem 12. II The observed 


value of T is a, and distribution of T does not depend on unknown parameters a and /3. 

















Therefore, T is a generalized pivotal variable for a, and can be used to construct a generalized 
confidence interval for a. 

Let 

2 

T* = T-a = rniy) -a. 

Then, T* is a generalized test variable for a, because i) the observed value of T* does not depend 
on any parameters, ii) the distribution function of T* is free from nuisance parameters and only 
depends on the parameter a, and iii) the distribution function of T* is an increasing function 
with respect to the parameter a, and so, the distribution of T* is stochastically decreasing in 
a. Therefore, the generalized p-value for testing the hypotheses in (12. 8 j) is given as 

p = P{T* <0\Ho) = P{T <ao), ( 2 . 11 ) 

and the generalized p-value for testing the hypotheses in (j2.9|) is given as 


p = 2 min {P (T > oq) ,P{T< ao)}- 


( 2 . 12 ) 


and 


The generalized confidence interval for a based on T, and the generalized p-values in (12.111) 
)2.12p can be computed using Monte Carlo simulation ()Weerahandill995l : 


Behboodian and .lafari 


20061 1 based on the following algorithm: 


Algorithm 2.1. For given ro,ri,... ,rn, 

1. Generate U ~ x^ 2 n) ^ ~ X( 2 n+ 2 )- 

2. Compute T in (j2.10p . 

3. Repeat steps 1 and 2 for a large number times, (say M = lOOOOj, and obtain the values 
Ti,..., Tm ■ 

4- Set Di = 1 if Ti < ao else Di = 0, I = 1 ,..., M. 

The 100(1 — 7 )% generalized confidence interval for a is ^(1-7/2)] ; where is 

the 5th percentile of Tfis. Also, the generalized p-value for testing the one-sided hypotheses in 
(j2.1ip is obtained by A- 


2.3 Real example 


Roberta (jl979l l gave monthly and annual maximal of one-hour mean concen tration of sulfur 


dioxide (in pphm) from Long Beach, California, for 1956 to 1974. IChanI (119981 ) showed that the 





























Weibull distribution is a reasonable model for this data set. 


Wu and Tsena ( 2006l l also study 


this data set. The upper record values for the month of October from the data are 


26, 27, 40, 41. 

The 95% confidence interval for the scale parameter a based on our generalized confidence 
interval with M = 10000 is obtained as (5.4869,39.9734). The 95% confidence interval for 
the shape parameter /3 in (j2.6l) is obtained as (0.6890,8.0462), and based on Wu and Tseng’s 
method in (12.71) is obtained as (0.6352, 7.7423). Also, the generalized p-value equals to 0.0227 
for testing the hypotheses in (12.8p with ao = 5. Therefore, the null hypothesis is rejected. 


Table 1: Empirical coverage probabilities and expected lengths of the generalized confidence 
interval for the parameter a with confidence level 0.95. 





13 


a 

n 

0.5 

1.0 

1.2 

1.5 

2.0 

3.0 

5.0 

Empirical 

1.0 

3 

0.951 

0.949 

0.953 

0.947 

0.947 

0.948 

0.948 

Coverage 


7 

0.952 

0.949 

0.950 

0.950 

0.951 

0.953 

0.952 



9 

0.951 

0.948 

0.953 

0.951 

0.948 

0.949 

0.950 



14 

0.945 

0.949 

0.950 

0.950 

0.954 

0.952 

0.952 


2.0 

3 

0.949 

0.952 

0.947 

0.949 

0.951 

0.950 

0.953 



7 

0.948 

0.953 

0.950 

0.946 

0.954 

0.948 

0.951 



9 

0.952 

0.948 

0.953 

0.950 

0.952 

0.953 

0.954 



14 

0.950 

0.946 

0.949 

0.951 

0.951 

0.952 

0.955 

Expected 

1.0 

3 

16.740 

3.581 

2.804 

2.155 

1.653 

1.211 

0.847 

Length 


7 

13.575 

3.198 

2.477 

1.942 

1.475 

1.041 

0.686 



9 

13.505 

3.138 

2.469 

1.918 

1.446 

1.008 

0.651 



14 

13.122 

3.082 

2.403 

1.854 

1.376 

0.943 

0.596 


2.0 

3 

33.516 

7.187 

5.579 

4.341 

3.323 

2.427 

1.704 



7 

27.960 

6.344 

4.999 

3.899 

2.960 

2.080 

1.364 



9 

27.342 

6.304 

4.935 

3.831 

2.890 

2.016 

1.302 



14 

26.626 

6.129 

4.779 

3.705 

2.757 

1.886 

1.191 


2.4 Simulation study 

We performed a simulation study in order to evaluate the accuracy of proposed methods for 
constructing confidence interval for the parameters of Weibull distribution. For this purpose, 
we generated n + 1 record values from a Weibull distribution, and considered a = 1,2. For 
the simulation with 10000 runs and different values of the shape parameter /3, the empirical 
coverage probabilities and expected lengths of the methods with the confidence coefficient 0.95 
were obtained. The results of our generalized confidence interval for inference on a using the 















algorithm 12.11 with M = 10000 are presented in Table [H and the results of our exact method 
(E) and the Wu method (W) for inference on /3 are given in Table [2l We can conclude that 

i. The empirical coverage probabilities of all methods are close to the confidence level 0.95. 

ii. The expected lengths of E and W increase when the parameter /3 increases. Additionally, 
the expected length of E is smaller than W especially when /3 is large. 

hi. The expected length of our generalized confidence interval for a decreases when the 
parameter (3 increases. Moreover, it is very large when /3 is small. 

iv. The expected lengths of all methods decrease when the number of records increases. 

V. The empirical coverage probabilities and expected lengths of W and E do not change 
when the parameter a changes. 


3 Joint confidence regions for the parameters 


Suppose Rq, i?i,..., Rn are the first n + 1 upper record values from a Weibull distribution with 
parameters a and /?. In this section, we presented a joint confidence region for the parameters 
a and (3. This is important because it can be used to find confidence bounds for any function of 
the parameters such as the reliability function R (t) = exp( —Eor more references abo ut 


the 

and 


pint confidence reg i on ba sed on records, reader can see 


Asgharzadeh et al.l ( 20111 ). 


Asgharzadeh and Abdil (j 2 nila 


u 


3.1 Asgharzadeh and Abdi method 


Asgharzadeh and Abdil (j2011bl ) present exact joint confidence regions for the 


Weibu ll distribution based on the record values using the idea presented by 


parameters of 


Wu and Tseng 


(1200a)- The following inequalities determine 100(1 — 7 )% joint confidence regions for a and 

/ 3 : 

logff^Ui + l) log((^)k2 + l 


Aj = 


log 


Rn 


Rn 


Rj-1 

2 


</ 3 < 


log 


Rn 

Rj — 1 


(3.1) 


,^(2n+2),(l+vT^)/2, 


< a < Rn 


X 


( 2 n+ 2 ),(l-vT^)/ 2 , 


for j = 1 ,..., n, where 


ki = F, 


k2 = F( 
























Table 2; Empirical coverage probabilities and expected lengths of the methods for constructing 


confidence interval for the parameter /? with confidence level 0.95 






/? 


a 

n 

Method 

0.5 

1.0 

1.2 

1.5 

2.0 

3.0 

5.0 

Empirical 

1.0 

3 

W 

0.950 

0.952 

0.952 

0.949 

0.949 

0.953 

0.947 

Coverage 



E 

0.949 

0.953 

0.953 

0.950 

0.948 

0.953 

0.946 



7 

W 

0.951 

0.949 

0.950 

0.948 

0.949 

0.949 

0.953 




E 

0.951 

0.950 

0.948 

0.950 

0.953 

0.953 

0.950 



9 

W 

0.949 

0.949 

0.945 

0.949 

0.947 

0.949 

0.950 




E 

0.948 

0.950 

0.948 

0.948 

0.949 

0.952 

0.949 



14 

W 

0.946 

0.949 

0.951 

0.951 

0.953 

0.949 

0.951 




E 

0.947 

0.948 

0.950 

0.951 

0.953 

0.952 

0.952 


2.0 

3 

W 

0.954 

0.955 

0.948 

0.950 

0.950 

0.952 

0.949 




E 

0.953 

0.953 

0.947 

0.949 

0.949 

0.952 

0.950 



7 

W 

0.950 

0.955 

0.947 

0.948 

0.952 

0.947 

0.948 




E 

0.949 

0.952 

0.951 

0.948 

0.952 

0.947 

0.950 



9 

W 

0.953 

0.948 

0.956 

0.950 

0.953 

0.951 

0.952 




E 

0.952 

0.948 

0.951 

0.951 

0.953 

0.951 

0.953 



14 

W 

0.948 

0.947 

0.949 

0.950 

0.951 

0.951 

0.952 




E 

0.950 

0.947 

0.949 

0.949 

0.953 

0.951 

0.955 

Expected 

1.0 

3 

W 

1.704 

3.431 

4.194 

5.279 

6.913 

10.396 

17.479 

Length 



E 

1.630 

3.285 

4.013 

5.041 

6.611 

9.937 

16.722 



7 

W 

0.932 

1.879 

2.224 

2.808 

3.752 

5.578 

9.276 




E 

0.853 

1.716 

2.038 

2.574 

3.437 

5.115 

8.499 



9 

W 

0.806 

1.603 

1.928 

2.412 

3.222 

4.797 

8.024 




E 

0.730 

1.450 

1.748 

2.185 

2.928 

4.352 

7.267 



14 

W 

0.625 

1.262 

1.509 

1.888 

2.505 

3.766 

6.263 




E 

0.558 

1.125 

1.343 

1.685 

2.236 

3.353 

5.590 


2.0 

3 

W 

1.713 

3.458 

4.156 

5.277 

6.856 

10.307 

16.998 




E 

1.638 

3.306 

3.967 

5.053 

6.560 

9.859 

16.266 



7 

W 

0.934 

1.866 

2.208 

2.822 

3.738 

5.629 

9.392 




E 

0.854 

1.710 

2.026 

2.589 

3.419 

5.151 

8.581 



9 

W 

0.808 

1.600 

1.923 

2.418 

3.189 

4.798 

8.032 




E 

0.733 

1.451 

1.743 

2.193 

2.890 

4.347 

7.276 




W 

0.628 

1.260 

1.496 

1.888 

2.523 

3.768 

6.302 




E 

0.560 

1.124 

1.338 

1.688 

2.250 

3.366 

5.624 




and is the 7 th percentile of the F distribution with a and b degrees of freedom. Note 


that for each j, we have a joint confidence region for a and /?. 


Asgharzadeh and Abdil () 2011 bl i 


found that in most cases A|^n+ij and A|^n+i_|_^j provide the smallest confidence areas, where [xj 
is the largest integer value smaller than x. 

3.2 A new joint confidence region 

From Theorem 12.11 V = 2/3^”^Qlog has a chi-square distribution with 2n degrees of 

freedom and V = 2(-^)^ has a chi-square distribution with 2n -|- 2 degrees of freedom, and U 
and V are independent. Therefore, an exact joint confidence region for the parameters a and 
j3 of Weibull distribution based on the record values can be given as 


B = 


2 2 
X(2n),(l-v7^)/2 ^ ^(2n),(l-HvT^)/2 


Rn 


R„ 


.^(2r^-h2),(l-h^/T=7)/2. 


2 Si=0 log R, 

< a < Rn 


(3.2) 


X 


(2n-H2),(l-^/T=7)/2, 


Remark 3.1. All record values are used in the proposed joint confidence region in but 

not in the proposed joint confidence regions in m- 

3.3 Real example 

Here, we consider the upper record values in the e xample given in Section 12.31 T herefore, the 


95% joint confidence regions for a and j3 based on 


Asgharzadeh and Abdil ( 2011bl i in (13.ip are 


Ai = |(a,/3) : 0.5826 </? < 11.9955, 41(0.1029)? < a < 41(1.1318)?} 

A 2 = {(a,/3) : 0.1646 < j3 < 6.4905, 41(0.1029)? < a < 41(1.1318)?} 

H 3 = |(a,/3) : 0.1720 </3 < 58.9824, 41(0.1029)? < a < 41(1.1318)?} 

and the 95% joint confidence region for a and fi in (|3.2p is 

R = |(a,/3) : 0.5305 </3 < 9.0277, 41(0.1029)? < a < 41(1.1318)?} . 


The plot of all joint confidence regions are given in Figure [TJ Also, the area of the joint 


confidence regions Ai, A 2 , A 3 , and B are 194.9723, 166.7113, 369.7654, and 172.5757, respec¬ 
tively. 
















90 rr 
80 - 
70 - 
60 - 



0 10 20 30 40 50 60 

P 



Figure 1: The plot of the joint conhdence regions Ai, A 2 , A 3 , and B. 


3.4 Simulation study 

We performed a similar simulation given in Section 


with considering a = 1, in order 


to compare the joint conhdence regions proposed bv lAseharzadeh and Abdil (1201 Ibl i and our 
joint conhdence region (B) in ()3.2I) . Here, we consider the conhdence areas Ay^r^^ and 


A 


n+i I ^ j because the coverage 


probabilities of all Ai ’s are close to the conhdence coefficient and 


Asgharzadeh and Abdil (j2nilbl ) found that in most cases these two conhdence areas provide 


the smallest conhdence areas. The empirical coverage probabilities and expected areas of the 
methods for the conhdence coefficient 95% are given in Table El We can conclude that 


1. The coverage probabilities of the all methods are close to the conhdence coefficient 0.95. 


2. The expect ed area of our method i s sma. 
methods bv lAsgharzadeh and Abdil (l2nilbl ) 


ler than the expected areas of the proposed 


3. The expected areas of all methods decrease when the number of records increases. 


4. The expected areas of all methods decrease when the parameter /3 increases. 

























Table 3: Empirical coverage probabilities of the methods for constructing joint confidence 


region for the parameters a and /3 with 7 = 0.05. 





/? 


n 

Region 

0.5 

1.0 

1.2 

1.5 

2.0 

3.0 

5.0 

Coverage 

4 

Ai 

0.950 

0.949 

0.949 

0.951 

0.954 

0.946 

0.950 

Probability 


A2 

0.949 

0.951 

0.950 

0.951 

0.953 

0.949 

0.952 



B 

0.949 

0.950 

0.949 

0.952 

0.954 

0.949 

0.950 


6 

Ai 


0.951 

0.952 

0.952 

0.954 

0.949 

0.949 



A2 


0.948 

0.950 

0.948 

0.953 

0.948 

0.952 



B 


0.949 

0.951 

0.951 

0.953 

0.950 

0.951 


9 

A2 


0.949 

0.950 

0.955 

0.949 

0.948 

0.953 



A3 


0.951 

0.951 

0.956 

0.949 

0.948 

0.949 



B 


0.952 

0.951 

0.953 

0.949 

0.948 

0.952 


14 

A3 

0.947 

0.950 

0.954 

0.948 

0.954 

0.951 

0.952 



A4 

0.946 

0.951 

0.952 

0.948 

0.951 

0.950 

0.952 



B 

0.948 

0.949 

0.952 

0.947 

0.951 

0.953 

0.952 


29 

Aq 


0.953 

0.950 

0.950 

0.948 

0.948 

0.950 



A7 


0.952 

0.951 

0.949 

0.948 

0.948 

0.950 



B 


0.955 

0.951 

0.953 

0.950 

0.948 

0.952 

Expected 

4 

Ai 

27.787 

8.548 

7.339 

6.331 

5.725 

5.330 

5.203 

Area 


A2 

30.020 

8.976 

7.651 

6.682 

5.989 

5.593 

5.504 



B 

22.985 

7.371 

6.388 

5.596 

5.099 

4.792 

4.713 


6 

Ai 

21.062 

6.036 

5.213 

4.576 

4.102 

3.783 

3.701 



A2 

20.035 

5.824 

5.046 

4.436 

3.985 

3.701 

3.648 



B 

14.551 

4.714 

4.144 

3.714 

3.399 

3.192 

3.162 


9 

A2 

14.631 

4.081 

3.533 

3.059 

2.756 

2.574 

2.510 



A3 

14.651 

4.086 

3.545 

3.077 

2.767 

2.585 

2.541 



B 

9.639 

3.137 

2.774 

2.471 

2.275 

2.160 

2.133 


14 

^3 

9.436 

2.702 

2.320 

2.035 

1.816 

1.701 

1.661 



Aa 

9.388 

2.686 

2.304 

2.035 

1.812 

1.707 

1.668 



B 

5.784 

1.999 

1.763 

1.599 

1.471 

1.405 

1.385 


29 

Aq 

4.244 

1.298 

1.124 

0.988 

0.905 

0.849 

0.828 



Aj 


1.291 

1.118 

0.985 

0.904 

0.850 

0.828 



B 


0.932 

0.838 

0.761 

0.719 

0.689 

0.678 


























Appendix. Generalized j)-value and generalized confidence in¬ 
terval 

Let X be a random variable whose distribution depends on a parameter of interest 9, and a 
nuisance parameter A. Let x denote the observed value of X. A generalized pivotal quantity 
for 0 is a random quantity denoted by T{X;x\0) that satisfies the following conditions: 

(i) The distribution of T{X;x;9) is free of any unknown parameters. 

(ii) The value of T{X;x;6) at X = x, i.e., T{x-,x;9) is free of the nuisance parameter A. 
Appropriate percentiles of T(X]x;6) form a confidence interval for 9. Specifically, if 

T{x;x;9) = 9, and denotes the 100(5 percentage point of T{X-,x]9) then (Ty/ 2 , Ti_^/ 2 ) 
is a 1 — 7 generalized confidence interval for 9. The percentiles can be found because, for a 
given X, the distribution of T(X; x] 9) does not depend on any unknown parameters. 

In the above setup, suppose we are interested in testing the hypotheses 

Ho:9<9o vs. Hi:9> 9o, (A.l) 

for a specified 0o- The generalized test variable, denoted by T*(X; x-,9), is defined as follows: 

(i) The value of T*(X; x;9) at X = x is free of any unknown parameters. 

(ii) The distribution of T*(X; x] 9) is stochastically monotone (i.e., stochastically increasing or 
stochastically decreasing) in 9 for any fixed x and A. 

(hi) The distribution of T*(X; x] 9) is free of any unknown parameters. 

Let t* = T* {x;x\9q)., the observed value of T*{X;x;9) at {X-,9) = {x-,9q). When the 
above conditions hold, the generalized p-value for testing the hypotheses in (lA.lh is defined as 

p = P{T*{X-,x-,9o)<r) (A. 2 ) 

where T*{X-,x-,9) is stochastically decreasing in 9. The test based on the generalized p-value 
rejects Hq when the generalized p-value is smaller than a nominal level 7 . However, the size 
and power of such a test may depend on the nuisance parameters. 
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